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Abstract 

This paper presents an empirical investigation of the intraday Brazilian stock market 
price fluctuations, considering q-Gaussian distributions that emerge from a non- 
extensive statistical mechanics. Our results show that, when returns are measured 
over intervals less than one hour, the empirical distributions are well fitted by q- 
Gaussians with exponential damped tails. Scaling behavior is also observed for these 
microscopic time intervals. We find that the time evolution of the distributions is 
according to a super diffusive q-Gaussian stationary process within a nonlinear 
Fokker-Planck equation. This regime breaks down due to the exponential fall-off of 
the tails, which in turn, governs the transient dynamics to the long-term macroscopic 
Gaussian regime. Our results suggest that this modeling provides a framework for 
the description of the dynamics of stock markets intraday price fluctuations. 
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1 Introduction 



The empirical probability distribution functions (PDFs) of price fluctuations 
of financial indices for different markets have been reported in the econophysics 
literature in recent years [1,2]. Although there has been a huge progress in the 
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statistical description of these fluctuations, a complete and consistent descrip- 
tion of the distributions and its dynamics in the high-frequency regime is still 
lacking. 



It is well known that high-frequency financial time series such as foreign ex- 
change rates or stock market returns are long-memory non-Gaussian processes. 
Such complex behavior calls for advanced theories in the field of statistical 
physics. Whence, we consider the recently proposed non-extensive statistical 
physics [3], a theory that generalizes the extensive classical statistical theory 
for strong dependent variables. A large number of studies [4,5,6,7] of financial 
markets have already employed the non-extensive statistics in their analysis. 



In this paper, we investigate the intraday dynamics of the BOVESPA index, 
the financial index of the Brazilian stock market, one of the largest emerging 
markets in the world. We model the distributions of the index price fluctu- 
ations by q-Gaussians, which are a class of stable distributions that emerge 
from the non-extensive approach, and where the parameter q measures the 
degree of the non-extensivity of the stochastic process. 



We find that for time horizons less than one hour, the probability distributions 
of price fluctuations are well fitted by q-Gaussians with damped exponential 
tails. A q-Gaussian scaling at the centre of the distributions with q*=1.75 is 
observed, which holds over these high frequency time scales due to the presence 
of strong correlations. The quasi-stationary q*-Gaussian regime breaks down 
due to the exponentially truncated tails and a new intraday correlated regime 
emerges at larger time scales, in which the q*-Gaussian central part of the 
distributions are consumed by the exponential tails. 



The dynamics in the short-time regime is in agreement with the prediction 
of a super diffusive q-Gaussian stationary process governed by a nonlinear 
Fokker-Planck equation (FPE) that models the correlated anomalous diffusion 
found for high frequency price fluctuations. Therefore, we present a coherent 
description encompassing non-extensive statistics and an evolution equation. 



This paper is organized as follows. In section 2, we describe the empirical 
observations for the high-frequency BOVESPA index. In section 3, we present 
the non-extensive statistical theory. In sections 4 and 5, these observations are 
analyzed by using the non-extensive approach. In section 6, we present some 
concluding remarks. 
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2 Empirical results 



Our results are based on analyzing BOVESPA index high-frequency 30-seconds 
interval data collected from November 2002 to June 2004, a period without 
major local market disturbances. This data consists of a set of 352500 prices 
(index values) p(t) which allows a fairly statistical analysis. 

From this database, we select a complete set of non-overlapping price fluctu- 
ations (returns): 

x T (t) = In p(t + t) - In p(t) (1) 
occurring in r-intraday interval. We define the normalized returns as: 

x?(t) = (x T (t)-fi T )/a T (2) 
where \i T and a T are respectively the mean and the standard deviation of x T (t). 

To characterize quantitatively the observed stochastic process, we measure the 
probability distribution P(x^) of index price fluctuations for r ranging from 
1 to 60 minutes. The number of data in each set decreases from the maximum 
value of 176250 (r = 1 min) to the minimum value of 29375 (r = 60 min). 

10% 
10 s - 

ioS 

I 10 3 - 



10 1 1 
10°- ■ 

-10 -5 5 10 

Normalized return 




Fig. 1. Binned non- normalized frequency distribution of Bovespa index r = 1 minute 
normalized returns in logarithmic scale (full squares) and comparison with the ex- 
pect distribution for a Gaussian process (open squares). 

In figure 1 we show the measured (non-normalized) probability distribution for 
t = 1 min normalized returns and compare with the expected distribution for 
the Gaussian process. The empirical distribution is clearly fat tailed, indicating 
that non-Gaussian time series analysis is required. 
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The scaling behavior of the distributions at coarser time scales has different 
regimes. At micro scales (typically shorter than one hour), correlations be- 
tween successive price changes are dominant due to informational flow between 
the financial trades. On the other hand, at macro scales (typically larger than 
one month) macroeconomic rules govern the market drift and correspond to 
a Gaussian regime. Mesoscopic time scales can be associated with a transient 
regime between the microscopic and macroscopic ones. 

It is still an open question the evolution of the PDFs at micro and meso 
time scales. In this paper, we investigate the scaling behavior of the return 
distributions at microscopic time scales. As shown in figure 2, the data collapse 
into the r = 1 min distribution, but this stationarity breaks down at r = 1 
hour. This means that the r = 1 min return distribution functional form is 
preserved for several micro time scales, therefore exhibiting a short-time non- 
Gaussian intraday quasi-stationary regime. 
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Fig. 2. Distributions of normalized returns of the Bovespa index sampled with differ- 
ent time scales (shown in the inset) (a) the data collapses into a single distribution, 
(b) the collapse breaks down at r=60 min time scale. 



We also investigate the presence of linear and non-linear dependence between 
successive r = 1 min normalized returns. In Figure 3(a) we show the linear 
autocorrelation R(At) =< x^ (t)x^ (t + At) > according to the time lag At 
of the measured returns. Employing the standard least-squares fit for time 
intervals r = 2-15 minutes, the semi-log plot indicates a short range exponen- 
tial behavior R(At) ~ exp(— At/T) with characteristic time scale T = 4 min. 
After At = 20 min the correlation is at the noise level. 
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Fig. 3. (a) semi-log plot of the autocorrelation function of r = 1 min normalized 
returns according to the temporal distance At of the measured returns. The solid 
line represents the linear fit from At = 2 min to 15 min. The autocorrelation reaches 
noise level after At = 20 min. (b) log-log plot of the covariance function of the 
amplitude of r = 1 min normalized returns according to temporal distance At of 
the measured returns. The solid line represents the linear fit from At = 2 min to 
60 min. 



In Figure 3(b) we show the covariance of the amplitude of the r = 1 min 
normalized returns C( At) =< \x? (t)\ \x?(t+At)\ > - < \x? (t)\ > 2 according 
to lag At. Employing the standard least-squares fit for time intervals r = 2- 
60 minutes, the log-log plot shows a power-law dependence C(At) ~ At~ v 
with 77 = 0.3 along all intraday time lags. This result signals a long- memory 
process, without a measurable temporal scale of decay of correlation within the 
analyzed regime. These findings are quantitatively similar to previous results 
for the U.S. market [8]. 



The anomalous diffusive process driven by the long and/or short range de- 
pendence of price fluctuations can be characterized by the scaling behavior 
o"(t) ~ r H of the standard deviation of the distributions. Through a log-log 
plot shown in figure 4 we estimate the scaling (Hurst) exponent H = 0.7. This 
non-trivial scaling exponent is the signature of a super diffusive behavior for 
the intraday time scales. Actually, it is a stylized fact that the intraday finan- 
cial market generates strong dependence in the consecutive price formation. 
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Fig. 4. Log-log plot of the standard deviation of the empirical return distributions 
as a function of time scale r. The solid line represents the linear fit from r = 1 min 
to r = 60 min. 

3 Non-extensive statistics 



The quasi-stationarity property of the microscopic time scale distributions 
shown in figure 2 suggests the search for stable distributions to model them. 
Sums of independent or weak dependent stochastic variables lead to Gaussian 
or to Levy limit regimes [9]. On the other hand, the presence of linear and 
non-linear dependence found for the one minute returns suggests the search 
for stationary distributions in a long memory environment. 

Recently, a non-extensive statistical approach that generalizes the standard 
Boltzmann-Gibbs theory has been developed, as an attempt to describe non- 
equilibrium regimes and/or systems with strong dependence. 

The cornerstone of this formalism is the entropy functional S q [p] = — 1 ii^M! dx 
where q < 3 is a real parameter representing the degree of non-extensivity of 
the functional and p(x) is a normalized probability density function. The q — > 
1 limit recovers the classical extensive property. The optimization of S q under 
natural constraints leads to q-Gaussian PDFs defined as: 



G q (x, t) = {l - (3(r) [(1 -q)(x- x(r)Y] } + ^ 



(3) 



where Z q {r) is a normalization constant, (3(t) is a scale parameter (fl^ 1 is 
proportional to the variance of the distribution) and the subindex + indicates 
that G q (x,r) = if the expression inside the brackets is non-positive. The 
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usual Gaussian distribution is recovered in the q — > 1 limit. 

An important property of the q-Gaussian (3) is that, with appropriate time- 
dependent parameters f3{r) and Z q (r), it is the invariant solution for a class 
of (correlated) anomalous diffusive processes governed by non-linear Fokker- 
Planck equation of the form [10]: 

= ~Y X \n*)P& r)\ + [P ^(x, r)] (4) 

where F(x) = a — b.x is a mean reverting force with 6 ^ 0. The explicit time 
dependence of (3(r) and Z q {f) are: 

/3(r)^ = P^^exp [-6(3 - q)(r - r )] + 

+2Db-\2 - q) [(3(ro)Z 2 q (r )} V {1 - exp [-6(3 - q)(r - r )]} (5) 

Z q (r) = Z q (r ) {(1 - A q )exp [-(r - r )(6(3 - q))} + AJ^ (6) 

where A g = 26 _1 (2 — g)D/?(r )Z^ 1 (ro). Equation (5) describes anomalous 
scaling of the inverse variance parameter (5 for q ^ 1. In the limit of weak 
reverting force and for time scales such that 6(3 — q) « (r — tq)^ 1 << 
2D( y 2 — q)(3 — q)f3(To)Z q ~ 1 (T ), one gets the result for the free diffusion process 
(F(x)=0) [11]: 

Super diffusion occurs for q > 1. This is the range of interest of the non- 
extensive parameter for application to our intraday data. Hence, parameter q 
represents the degree of the resulting anomalous diffusion from the underlying 
interaction among financial trades. 

We model the normalized returns at time scale r by the stochastic variable 
x in Eq. (4). Thus, the solutions G q (x, r) describe the empirical distribution 
P«). 

The mean-reverting parameter a controls the average equilibrium value of 
the stochastic variable and does not affect the diffusion properties. In our 
application for normalized returns, a is set equal to zero. On the other hand, 
from (5), the mean- reverting rate parameter 6 is crucial for the attainment of 
a long-time equilibrium solution, but it is worth noting that, as b gets smaller, 
the spread of the distributions holds over longer periods and the equilibrium 
solution may not be observable. 

Substituting the explicit expression of G q (x, r) in (4) it can be shown [12] 
that the stochastic variables, whose addition has the q-Gaussian as a limit 
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distribution, are characterized by non-null linear correlation (for b ^ 0) and 
non-null higher-order correlations (for q ^ 1). In our application, these step 
variables are the one minute normalized returns, which, according to the em- 
pirical findings, are consistent with this modeling with g ^ 1 and b ^ 0, for 
microscopic time scales. 



4 Non-extensive modeling of the intraday return distributions 

It is apparent from the results sketched in the previous section that the q- 
Gaussians are good candidates for describing the anomalous diffusion of fi- 
nancial market indices. Since the correlations are stronger for shorter trading 
time intervals, we focus on the microscopic time scales less than one hour (1 
min < t < 60 min). 
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Fig. 5. (a) Probability distribution function of normalized returns for time interval 
r = 1 min: BOVESPA index data (black squares) and best fit of the data with Eq 
(3): (a) range xf G [-10, + 10]; best fit q=1.64 and /?=3.36. (b) range G [-5, +5]; 
best fit q*=1.75 and /3=4.47. 

We determine the parameters of the model by minimizing the mean square 
deviation between the empirical distribution P(x^) and the q-Gaussian dis- 
tribution form G q (x,r). This method was applied throughout this paper. The 
free parameters are q and (5 controlling respectively the shape and the width 
of the q-Gaussian. In Figure 5(a) we plot the binned empirical distribution of 
normalized r = 1 min returns and the optimal fit. We find that the empiri- 
cal distribution is well fitted with a q-Gaussian curve within the range Xi E 
[-5, +5] and it is overestimated outside this range. The optimal non-extensive 
parameters are q = 1.64 and (3 = 3.36of , where <j\ = 5.7. 10~ 4 . 
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We now search for the optimal fit only for the bulk of the r = 1 min returns, 
defined as returns with magnitude less than or equal 5a±. Figure 5(b) shows 
that a good agreement with the q-Gaussian profile is obtained when ^^(r)! < 
5 with q* = 1.75 and (3 = 4.47/(7?. 

Comparing the above results, we find that the inclusion of the tails of the 
empirical distributions in the fit has the effect of diminish the optimal q- 
value. Nevertheless, the fit at the tails are not satisfactory. Then, for our 
non-extensive modeling, we consider only the bulk of the distribution. 

To analyze the stability of the market through the elapse time of observation 
of our data, we estimate the local parameter q and investigate if it is highly 
fluctuating. In figure 6, we show the time series of the local q parameters 
considering the bulk of r = 1 min returns for each month. From this figure, 
they are roughly constant over the months and are consistent with the optimal 
value q*=1.75 for the entire data set. Hence, it is possible to assume that this 
parameter is an invariant of the market characterizing the micro time scales. 
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Fig. 6. Time series of optimal parameter q for the bulk of r = 1 min returns collected 
in each month of the data set. The horizontal line represents the optimal parameter 
q*=1.75 for the whole period of the analyzed data. 

In figure 5(a), it was shown a clear deviation of the tails of the distribution 
from the q-Gaussian profile. To investigate the tails, cumulative distributions 
P cum (x) of the positive and negative tails for normalized r = 1 min returns 
are plotted in a semi-log space, as shown in figure 7. The linear dependence of 
the graph suggests that the r = 1 min returns exhibit an exponential decay 
in the tails. The straight lines are almost parallel to each other, indicating 
that both tails follow approximately the same exponential law P CU m(%) ~ 
exp(— (\x\ — Xc)H) with £ = Aa x . 
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Fig. 7. Semi-log plot of the cumulative distribution of positive and negative tails for 
r = 1 min normalized return distribution. 

This asymptotic exponential behavior for large returns also emerges as a ro- 
bust characteristic of the intraday distributions at larger time scales, as shown 
in figure 8, where the distributions exhibit a tent-shape form, when plotted in 
a semi- log space. 
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Fig. 8. Semi-log plot of the probability distribution of Bovespa index normalized 
returns in the range € [-5, +5]. For clarity, the distributions for successive r are 
shifted up by a factor of 10 each. The inset shows the time scales. 

The q-Gaussians have sharp peaks and power-law tails. For q > 5/3, the 
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q-Gaussian has an infinite second moment. However, for large returns, the 
empirical distribution crosses over to simple exponential decay. The exponen- 
tial fall-off implies that the second moment is finite. These results suggest that 
the distributions of high-frequency returns can be modeled by exponentially 
truncated q-Gaussian of the form: 

Pq( x ) = CP q (x), para \x\ < xq 

(8) 

Pq(x) = CP q (x)exp (- l a l~ ac ) , para \x\ > x c 

In Figure 9 we show that Eq. (8) very well reproduces the distribution of r 
= 1 min normalized returns of the Bovespa index with the fitting parameters 
q* = 1.75,(3 = 4.47/of,£ = 4<7i and x c = 4<7i. Note that the value of x c is 
consistent with the region \x±\ > 5a"! considered for the tails. 
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Fig. 9. Semi-log plot of the r = 1 min normalized probability distribution (open 
squares) and the best-fit using an exponentially truncated q-Gaussian given by Eq. 
(8) with parameters q* = 1.75, /J=4.47, £ = 4 and x c = 4. 

In section 2 it was shown through Figures 3-5 the presence of significant corre- 
lations of returns at microscopic time lags. It is also important to test the null 
hypothesis that the observed correlations cause the persistent fat tails of the 
high frequency distributions. In that case, destroying all correlations would 
almost vanish the fat tails for larger time scales. With this aim, we random- 
ized the empirical series of r = 1 min returns, by shuffling them, creating an 
artificial series. In table 1, we show a comparison between the effective best-fit 
parameter q for the empirical (EMP) distributions and for distributions of the 
T-scale returns generated by the aggregation of the r = 1 min shuffled returns 
(SHU). While the randomized data have reached the Gaussian regime (q=l) 
at microscopic time scales, this regime is far from been reached by the real 
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data in the intraday time scales. This result supports that the scaling prop- 
erty observed in figure 2(a) and the associated slow convergence to the stable 
Gaussian distribution are caused by the nontrivial correlations in the market 
data. 



Time scale r (min) 


EMP 


SHU 


1 


1.75 


1.75 


2 


1.75 


1.65 


4 


1.73 


1.50 


8 


1.64 


1.33 


16 


1.52 


1.24 


30 


1.51 


1.18 


45 


1.49 


1.07 


60 


1.48 


1.06 



Table 1 

Effective q-parameters for return distributions at several microscopic time scales, 
for the original empirical r = 1 min returns (EMP) and for the artificial shuffled 
r = 1 min returns (SHU). The effective parameter q was adjusted for \x T \ < 5a T , 
which includes partially the tails of the distributions when r > 1 min. 



5 Time evolution of the intraday return distributions 

In this section, we investigate the stationarity of the central part of the return 
distributions at microscopic time scales, shown in figure 2(a). We now con- 
sider the empirical distributions of rescaled returns, this is, detrended original 
returns in units of <j\. 

x?(t) = (x T (t) - ix T )lo x (9) 

We consider the r = 1 min optimal non-extensive parameter q*=1.75 and 
obtain the best-fit parameter (3{t) for the central part of the rescaled return 
distributions at coarser time scales. In figure 10, we plot the distributions of 
figure 8 in rescaled variables. It is shown that the optimal q*-Gaussians well re- 
produce the bulk of the r-min rescaled return distributions in the microscopic 
regime. 

The results presented so far lead to the hypothesis of a quasi-stationary q- 
Gaussian regime for microscopic time horizon. The presence of the exponential 
tails do not affect the q-Gaussian scaling at the centre of the distributions, 
which holds over the microscopic time scales due to the presence of strong 
correlations. 
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Fig. 10. Rescaled return distribution for microscopic time scales and optimal 
q*-Gaussian fit with parameter q* =1.75 : (a) r = 1, 2, 3 min (b) r = 5, 10, 
15 min. Distributions are shifted up by a factor 10 each, for best visualization. 

We now investigate the super diffusive time scaling behavior of the empirical 
distributions along the same lines of [13] . We test the observed behavior with 
the time evolution predicted by the non-linear Fokker Planck Eq. (4) with 
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Fig. 11. Time evolution of parameter (5 as a function of time scale r: Bovespa index 
data (black squares) and best fit (line) using Eq. (5) with q*=1.75 and tq= 1 min. 
The estimated model parameters are D=0.0286 ± 0.002 and b= 0.025 ± 0.003. 

In figure 11, we model f3(r) with Eq. (5), where q=q* and r = 1 min. There 
are two free coefficients left, namely the diffusion constant D and the mean- 
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reversion rate b. We show that the empirical values can be well-fitted by the 
theoretical values, for particular estimated coefficients D = 0.0286 ±0.002 and 
b = 0.025 ± 0.003. These numerical values are similar to previous results for 
the U.S. market [13]. Therefore, our results for the intraday Brazilian market 
indicate a time evolution according to the nonlinear Fokker-Planck form. 



6 Discussion 

The data collapse shown in figure 2(a) reveals the existence of an invariant 
time scaling master curve for the description of high-frequency returns. The 
observation of invariant quantities in the economic systems, as in the physical 
systems, denotes the existence of conservation laws or same rules governing 
these processes. 

Moreover, it is known that there is strong time dependence among the tick-by- 
tick price formation by the financial market traders. In this work, we model 
the high-frequency returns of BOVESPA index by considering the q-Gaussian 
distributions, which are stationary solutions for a class of anomalous diffusive 
processes driven by strong correlations. 

Adjusting only the central part of the one minute distribution with the q- 
Gaussian distribution function, we obtain the optimal non-extensive param- 
eter q*=1.75. This value is remarkably constant over the two-year period of 
our data, as shown in figure 6. Similar values (q=1.77 and q=2) were obtained 
for the central part of the return distributions for the Japanese [6] and Ko- 
rean [14] stock indices, respectively. On the other hand, this value is bigger 
than previous ones reported in the literature for the U.S. market [15], which 
were obtained adjusting the whole range of the empirical returns. Actually, as 
shown in Figure 5, the inclusion of the damped tails of the distributions in the 
fit has the effect of diminish the optimal q value, but the observed exponential 
fall-off behavior of our empirical data disable the modeling with q-Gaussians 
at the tails, due to their power-law character. 

Figure 9 shows that Eq. (8) very well reproduces the distribution of r = 
1 min normalized returns of the BOVESPA index. Exponential tails in the 
foreign exchange and stock return distributions have been reported for some 
world markets [16,17,18,19,20] for micro and meso time scales, although several 
authors have reported power-law decays [21,22,23], including those employing 
the non-extensive approach [4,5,7]. Our empirical distributions consist of a 
significant number of data which allows a fairly statistical analysis including 
the tails. The exponential tails may due to finite-size effects of the financial 
market, which prevents the development of scaling invariant power-law tails. 
This effect is expected to be stronger in the case of emerging markets. This 
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result is also in according to previous findings for the daily return distribution 
tail of the BOVESPA index [18]. 

We have shown that the index time series in the high frequency regime is a 
long memory process where the Hurst exponent is significantly greater than 
1/2. This long-memory leads to stationary leptokurtic distributions with the 
observed scaling behavior of /3(t) for microscopic time scales. The fast decay 
of P(t) signals this super- diffusive character of the market dynamics. 

The anomalous scaling (5) with two free parameters, D and b, reasonably fits 
the time evolution of the PDFs at micro scales. This modeling is able to predict 
the super diffusive behavior of this quasi- stationary regime. The estimated 
value for parameter b is small, but significantly different from zero. This result 
is consistent with non-null linear correlation of normalized returns for time lags 
At < 20 minutes shown in figure 3(a) and a lifetime of the regime towards the 
limit distribution that encompass the microscopic time scales. Recent analysis 
of the S&P500 signal [5] also shown that the first Kramers-Moyal coefficient 
of the FPE is almost zero, indicating almost no restoring force. The stable 
non-extensive parameter q* > 1 in Eq. (4) implies the existence of higher- 
order correlations [12], modeling the non-linear memory effects illustrated in 
figure 3(b). Hence, we conclude that the non-linear Fokker-Planck equation (4) 
captures the main features of the dynamical evolution of stock price returns 
at this time horizon. 

Another dynamical foundation of non-extensive statistics have been applied to 
financial indices [4,5], in which the price fluctuations are described by Brown- 
ian motion subject to a power-law restoring force and a Gaussian white noise 
with a slow varying amplitude. This leads to linear FPE for temporal evo- 
lution of the distribution function with time-dependent diffusive coefficient. 
However, if we take into account the strong correlations of the price changes at 
microscopic time scales, the description by means of a non-linear FPE seems 
to be more appropriate. 

The q-Gaussian regime breaks down due to the exponentially truncated tails 
that prevent the attainment of the equilibrium solution. On the other hand, 
the long-range correlations exhibited by the intraday diffusive process further 
delay the convergence to the long time Gaussian regime, giving rise to a new 
intermediate regime at mesoscopic time scales, in which the q-Gaussian central 
part of the distributions are consumed by the exponential tails. This result is 
in accordance with previous findings [17] for probability distributions of stock 
returns at mesoscopic time horizons that follow an exponential function. 

In conclusion, we found that the intraday stock price fluctuations in Brazil 
have an intermediate behavior between the non-extensive q-Gaussian and the 
extensive Gaussian regimes characterized by r-dependent exponential tails 
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of the return distributions. We have shown that a model based on a non- 
extensive approach that captures the strong dependence of the stochastic 
variable robustly accounts for the observed time evolution of the return dis- 
tributions at micro time horizons and reproduces the crossover between micro 
and meso time lags. This exponentially damped, non-extensive, dynamical be- 
havior should provide a framework to investigate other high-frequency stock 
market time series. 
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